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In the previous paper on this tòpic it was shown how, for a pulse of arbitrary shape and duration, 
the drive frequency can be analytically optimized to maximize the amplitude of the population 
oscillations between the selected two levels in a general many level quantum system. It was shown 
how the Standard Rabi theory can be extended beyond the simple two-level systems. Now, in order to 
achieve the quickest and as complete as possible population transfer between two pre-selected levels, 
driving pulse should be tailored so that it produces only a single half-oscillation of the population. 

, In this paper, this second (and final) step towards the controlled population transfer using modified 

^"") ■ (i.e. many level system) Rabi oscillations is discussed. The results presented herein can be regarded 

as an extension of the Standard 7r-pulse theory - also strictly vàlid only in the two level systems - to 

, the coherently driven population oscillations in general many level systems. 
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I. INTRODUCTION 



During the past 20 years a number of methods has been devised for state selective preparation and manipulation of 
discrete- level quantum systems [l|, |2j, y, |J, |5j . However, simple population oscillations, induced by a resonant driving 
pulse have received negligible attention as a prospective population manipulation method. This might be attributed 
• to two reasons. The first is that Rabi theory is based on the rotating wave approximation (RWA), and all attempts 
to generalize it without RWA (e.g. @, 0,0) are mathematically very involved. The second is that no attempt has 
i been made to analytically generalize the original Rabi theory beyond the two-level systems. 

In this paper an analytic extension of Rabi theory to transitions in many-level systems is presented. The aim is to 
'design' a driving pulse of the form: 

i ; F(t) = F m(t) cos(w(í) í) (1) 

c ■ 

by establishing analytical optimization relations between its parameters: maximum pulse amplitude Fq, pulse envelope 
shape m(t), and time dependent carrier frequency ui(t). The goal of this enterprize is twofold: the first is to achieve 
as complete as possible transfer of population between two selected states of the system; the second is to make this 
transfer as ràpid as possible. These two requirements, however, are conflicting: population transfer can be accelerated 
by using a more intense drive, but at the same time a stronger drive increases involvement of remaining system levels 
p\ ' in population dynamics hence deteriorating population transfer between a selected pair of levels. 

In the previous paper on this tòpic it was shown how, for a pulse of arbitrary shape and duration, the drive 
frequency can be analytically optimized to maximize the amplitude of the population oscillations between the selected 
two levels in a general many level quantum system. It was shown how the Standard Rabi theory can be extended 
beyond the simple two-level systems. Now, in order to achieve the quickest and as complete as possible population 
transfer between two pre-selected levels, driving pulse should be tailored so that it produces only a single half- 
oscillation of the population. In this paper, this second (and final) step towards the controlled population transfer 
using modified (i.e. many level system) Rabi oscillations is discussed. The results presented herein can be regarded 
as an extension of the Standard 7r-pulse theory (see e.g. ^(j) " a l so strictly vàlid only in the two level systems - to the 
coherently driven population oscillations in general many level systems. 



II. THEORETICAL ANALYSIS 



All the calculations in this section are done in a system of units in which h = 1. 
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A. Calculation setup 

A quantum system with N discrete stationary levels with energies Ei (i = 1, ...,N) is considcrcd. The system is 
driven by a time dependent perturbation given in Eq. l(T|). In the interaction picture, the dynamics of the system 
obeys the Schroedinger equation: 

|a(í) = -iV(í)a(í), (2) 

where a(t) is a vector of time-dependent expansion coefhcients oi(í), ajv(í)- The NxN Matrix V(í) describes 
interaction between the system and perturbation. Explicitly, its elements are given by: 

V i0 {t) = _2^M. m ( í )( e ^·("(*)-^i)< + é-itjMQ+wa) ty ^ 

Hij is transition moment between the i-th and the j-th levels induced by the perturbation. s,y = sign(Ei — Ej) and 
Lüij = \Ei — Ej\ are respectively the sign and the magnitude of the resonant frequency for the transition between the 
i-th and the j-th level. 

The aim is to induce population transfer between two arbitrarily selected levels, designated by a and (3, directly 
coupled by the perturbation (i.e. such that fi a p 0). To simplify equations, the time variable t is re-scaled to r, with 
transformation between the two given by: 

dr = *2tE0 m (t)dt. (4) 



Then with following substitutions: 



Eq. (0) transforms into: 



where: 



/•íj(t) = s i:j — (w(í)-Wy) (5) 

2 

g ió {r) = s i:j - (w(í)+Wy) (6) 

*W - ^SÉt(r) (7) 

Rii - ^ (8) 

Ma/3 



A a(r) = -iW(T)a(r), (9) 



Wij{r) = R ij (e ifi >l T)xlT) + e-^^W). (10) 



Initial conditions for the problem of selective population transfer comprise complete population initially (at t = 
t = 0) contained in only one of the selected levels, either a or (3. The other selected level, as well as all the remaining 
N-2 'perturbing' levels of the system are unpopulated at this time. 

Population evolution IL,(t) of the i-th level is determined from Ilj(t) = |ai(í)| 2 . 



B. Rabi-like population transfer in a three level system 

It was demonstrated in the previous paper on this tòpic that the analytical extension of the Rabi oscillations 
thcory bcyond two-level systems is anchored in the analysis of the simplest of the 'many-level' systems - a three level 
one. Hence, in this section the impact of the single additional level on the population transfer period is discussed: 
beyond the 'selected' levels a and /?, the system now contains one additional 'perturbing' level, designated with index 
p. The only requirements on the system infernal structure are that ^apil^pp ^ and [i ap — 0. While the first two 
requirements are necessary, the last one does not reduce the generality of the final results to any significant extent 
and is introduced for calculational convenience exclusively. 
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For the observed three-level system, the dynamical equation (JÜJ) reduces to: 



d_ 

dr 



a a (r) 

a p{ T ) 
a p (r) 



(11) 



ç'/o/J^Wt) _|_ e -ig a /3(r)x(T) g 

e — if a p(T)x(r) _|_ e ig a p(T)x(r) q Rp p (e~ ií ^ í - T ^ xí - T ^ + e 19 ?^^^) 













a a (r) 








_a p (r)_ 



1. Recapitulation: minimizing the impact of the perturbing level 

As it was shown in Q, Eq. when the driving frequency is near the resonant value for the transition a «-> f3, 

the following expression can be obtained for the population dynamics of level p: 



a p (T)*-a p (T)(* 0P i _ ) 



where 



&/3p 



- A /3 P 0)' 



(12) 



U(3 p — Lü a p 



(13) 



Put in words, with the conditions mentioned, the dynamics of the level p parametrically depends on the dynamics 
of the level (3 to which it is coupled. The relation between the amplitudes of the population oscillations for levels p 
and [3 follows directly from the above expression: 



Up « e /3 p(T)n /3 (T) 



where: 



e l3p 



A/3 P (r). 



(14) 



(15) 



Note that, as A ( g p (r) is generally very small and |to(t) < 1, that parameter <jp p actually determines the effective 
strength of applied perturbation: if a^ p << 1, then dynamical impact of level p is negligible and perturbation may 
be considered weak; if a 2 ^ p ~ 1, perturbation is very strong. 

Further, the requirement of the minimization of the dynamical impact of the perturbing level p on the transition 
a <-> (3 leads to the following equation for the optimized dynamics of the (a,(3) subsystem: 

(16) 

where the two- level state vector (b a (r), bp (r)) is merely the unitary transformed vector of the subsystem (a,/3): 

(17) 



d 




= —i 


1 






d7 






1 







"6a (t)" 


= e -iÀ(r) 


a Q (r) 









Note that the precise form of the real transformation matrix A(r) is irrelevant here as it has no impact on the 
population dynamics. 
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The optimization procedure produces the analytic expression for the chirp of the driving frequency, which in the 
lowest order of approximation (suitable for all but the most intensive perturbations) amovmts: 

u)(t) = ui a a + (s8 a sp p )(ui8p- uiap) — — ai\ í (m(t')Y 'd£ (18) 

W Q /3 + U!/3p yy t J 

and from the Eq l(ltj|) it is found that the time 6 required for a single population transfer between levels a and (3, 
determined from the fundamental relation of the 7r-pulse theory: 

jf*-f (19) 

equals (in units of r): 

= | (20) 

As was shown in the 'exact' numerical solution to the Eq. <| X X|) indeed maximizes the population oscillations 
in the a — (3 subsystem. However, as will be discussed below, the predicted value for the period of the population 
oscillations (Eq. 1)20(1) is somewhat smaller than the correct one, with the discrepancy increasing with the increasing 
population leak to the level p. In the following section this issue is resolved and the corrected analytical expression 
for determination of the population transfer (or oscillation) period is obtained. 

2. Patching the population conservation of the total systera 

To start the following discussion, notice that the optimized solution for the population transfer between levels a 
and (3 (described by the Eq. ((16(1 ) is unfortunately too good to be true. Namely, its serious drawback lays in the fact 
that the leak of the population from the (ct,/3) subsystem into the perturbing level p goes by completely unnoticed! 

Formally, the root of the problem hidcs in the fact that the mathematical trick which enabled decoupling of the 
level p dynamics from the rest of the system (i.e. the step between Eq.(18) and Eq.(19) in JjJ) destroys the unitarity 
of the full dynamical equation for the three level system, Eq. . The consequence is that the dynamical equation 
for the u — (3 subsystem, Eq. 1)16(1 itself claims to be unitary, keeping the population of that subsystem fully conserved. 
This is clearly impossible, as level p does indeed capture some population - the exact amount given by Eci. ((14|) . 

This malfunction caused by the decoupling procedure unfortunately cannot be remedied within the decoupling 
procedure itself - the patch has to be provided by an independent approach. To do this, the following argument is 
used: since it is the equation for the dynamics of level (3 that changes due to the decoupling procedure and consequently 
causes the breakdown of the population conservation, it is only the dynamical equation for the level (3 that has to be 
modified; then, as the dynamics of the level (3 is governed by the elements in the lower row of the dynamical matrix 
in Eq. ((16(l . a particular ansatz intervention precisely into these elements might help rectify the overall population 
dynamics of the whole three level system. Hence, the correction is sought in the following form: 



d 




= — i 


1 






~dr 


.6/3 (r). 




C(t) *£(t)_ 




bp{r)_ 



(21) 



where Ç\t\ and £[t] are real non-negative functions. Such an ansatz does not interfere with the phase-fitting 
effect of the driving frequency optimization and Eq. (|T%)| - forged by the decoupling procedure - which maximizes 
the population oscillations in the a — (3 subsystem, is left unharmed. Instead, it merely enables phase-independent 
modification of the amplitudes of a and (3 populations. 

Expressing requirement of population conservation in the total three-level system as: 

d\b a (r)\ 2 + d\b ( T )\ 2 + d\b p (r)\ 2 = , (22) 

splitting the phase and amplitude contributions of the three wave function projections on the three stationary 
states: 
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b a (r) = BaMe^-M , 

bf}(T) = B p {r)é*^ , (23) 
b p (r) = B p (T)e^) } 

and using the known relation between the populations of levels (3 and p, Ea. (|14|l . the following result is obtained: 



B a (r) dB a + (1 + e / , p (r))S /3 (r) dfy + -^(r) 2 de Pp = 



(24) 



Now the corrected equation for the a — (3 subsystem, Eq. (|21|) , can be used to eliminate dB a and dBp and introduce 
C[r] and £[t] in their stead: 



dB n 



-iBp 7m(e í( ^ (r) -^ (T)) ) dr , 
dB^a = -í(c(t) B (t) /m(e- l( ^ (T) -*° (T)) ) + £(r) 5,3 (^^(e^^)) dr 

Finally, taking together Ea. f2*4^ and Ea. lf2*5|l it is found that: 



(25) 



B a (r) (C(r) - 1 + £ ^ (t) ) sm(^(r) - a (r)) + B fi {r) ($(r) + ^ ln (l + ^Wj'Jsm^W) = 0. (26) 
Since for population oscillations Bp and -B Q are 180° out of phase, this condition can be satished only if: 

1 



C(r) = 



1 + £/3 P (r) 



(27) 



3. Population oscillation period modified 

Hcnce, the correct dynamical equation for the a — (3 subsystem which both maximizes the population oscillation 
amplitudes of these two levels as well as proper ly conserves the overall population of the three-level (a — (3 — p) system 
is: 



d 




= —i 


dï 


b (r)_ 





b a (r) 



(28) 



A simple extension of this result to the general many-level system (in which level a also has some perturbing levels 
jointly designated by q - attached to it) yields the total corrected dynamical equation for such a system: 



d 

dï 



b a (T) 



i-±\n{l + e aq (T)Y 



l+e Qg (r) 



b a {r) 



(29) 



Now to finalize the calculation of the corrected population transfer period the following procedure is administered. 
First, the time variable is transformed r — > ip according to: 



dr = n((p)d(p 



(30) 



The goal of this variable transformation is to produce, in the new time variable (/?, the closed dynamical equations 
for b a (ip) and bp(ip) describing the dynamics which is as close as possible to the simple harmònic oscillation. Second, 
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and to that end, the transformation Ea. H30(l is introduced into Ea . (|29|) . the resulting relation is differentiated with 
respect to ip and all but the lowcst order terms in the small parameters e aq (r) and e.p p (r) are kept. Hence the following 
result is established: 

d 2 b a {ip) + 1 d A (l + e aq {íp)) 3 (l + ep p {tp)) \ db a (ip) + k? & = Q 

dtp 2 2dip\ k 2 /dtp (1 + e aq ((p)) (l + epp(tpj) 

d 2 b/3(tp) , 1 d / (l + e ag (y))(l + e/3p(y)) 3 \ dbpjtp) | tf_ bp(ip)=Q (31) 



dtp 2 2dip\ k 2 ) dtp (l + £aq(?p)) (l + ep P (<p)) 

In the third and final step, the appropriate value of the free parameter n{tp) is selected 



K ^) 2 = (l + e aq (tp))(l + e Pp &)) (32) 

which transforms 

^ + ^(Mi + ^(,)))^) + *.M-o 

^^M 1+ *W»« + W^» (33) 

Both these equations are similar to the damped oscillator equation. For negligible damping (e aq (ip), ep p (tp) << 1), 
they reduce to the harmònic oscillator equations, in which case the population transfer time of tt/2 is obtained in 
the variable (p. In the the original time coordinate, r, the corrected time 6 required for a single population transfer 
between the levels a and (3 is then obtained from: 



Jo 



dT n (34) 



v/(l + e aï (r))(l + e^(r)) 2 



Note the difference between this result, and the result in Eq. (|19fl obtained from the Standard 7r-pulse theory: in 
the lowest order approximation, the corrected population transfer time is shorter then the one obtained from Eq. 119|) 
by an order of ^(e QÇ (r) + epp(r)). 

On the other hand, taking into consideration the damping factor in the Ea. (|33|l . approximating: 



— ln(l + e aq ^ p ((p)) w —e aq ^ p (ip>) (35) 

and using the damped oscillator theory |llj|fp.246). an increase in the damping (expressed as an increase in t aq (r) 
and e^p(r)) leads to the an increase in the population time transfer by an order of \-§^^ 1 aq p p { T )- As this correction 
is an order of magnitude smaller than the correction obtained from Ea. H34|l . it can be neglected for all practical 
purposcs. Furthermore, the additional corrections to the transfer period due to the neglected higher order elements in 
the Eq.J2U are of the same order of magnitude as this correction due to the first derivative component, and as they are 
impossible to obtain analytically, the value of this whole second order correction for the case of strong perturbations 
is somewhat shaky. However, this is not an issue, as the whole optimization theory developed in 9] - and on whose 
applicability the results of the above analysis hinge - assumes rather modest perturbations, and is not even expected 
to work properly for the extreme vàlues. 



III. NUMERICAL SIMULATIONS 



In this section, numerical simulations of systcm dynamics for unoptimized and fully optimized driving pulsc of the 
form of Eq. JQ) are presented and compared. Here, unoptimized driving pulse is the one with driving frequency equal 
to the pure resonant frequency between the two levels selected for population transfer (u>(i) = to a p) and with pulse 
duration T determined according to the Standard 7r-pulse theory relation, Eq. (|19|l . On the other hand, the parameters 
of the fully optimized pulse are determined from Eq. i|18fl and Ea. H34|) . 
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FIG. 1: In all of the numerical examples the same pulse form was used - m(t) = sin(í2í) 2 , but with different vàlues of maximum 
intensity parameter (Fo) and different total pulse duration, T. The respective vàlues are quoted in each particular example. 

In all cases, a three-level system is considered, with the following system parameters (a.u. = atòmic units): 
upa — 0.017671 a.u., sp a — 1, jjLp a = 0.073 a.u.; ujp p = 0.017611 a.u., Sf} p = — 1, jjL@ p = 0.098 a.u.. These 
system parameters correspond to the three ro-vibrational levels of the HF molecule in the ground electrònic state: 
a = (v — 0,j = 2,m = 0), /3 = (v = l,j = 1, m = 0), p = (v = 2,j = 2,m = 0). The pulse shape in all of the 
examples is m(t) = sin(ílí) 2 as shown in Fig 1. 

A. Population oscillations 

As was demonstrated in @ , frequency optimization minimizes the impact of the perturbing levels on the amplitude 
of the population oscillations. In this subsection, the necessity of the inclusion of additional correction Ea . (|34|1 for the 
population transfer time - alongside the correction for the driving frequency - will be demonstrated. Also, the validity 
and the limitations of the analytically obtained expression for this correction will be discussed. 

1. Legitimate perturbation 

Fig. 2 presents the dynamics of the system subjected to the external drive of limiting intensity, cr| p = 0.05, corre- 
sponding to the Fo = 2.80534* 10~ 4 a.u.. It is just strong enough to noticeably (albeit not significantly) distort the 
pure resonant oscillations, but at the same time weak enough so that the theory developed in [í| and further in this 
paper provides the full and precise quantitative corrections. 

Pulse duration determined according to the Standard 7r-pulse theory expression, Eq.^üJ) is T v = 3077832a. u., 
whereas the optimized one, obtained from Ea. (|34|l is T opt = 3126029a. u.. The pulse is aimed at producing five 
complete population oscillations. 

Three cases of dynamics arc prcscntcd: Fig. 2.a shows the unoptimizcd dynamics; Fig 2.b shows the 'semi-optimizcd' 
dynamics, with optimized driving frequency, but unoptimized population transfer period; finally, Fig. 2.c shows the 
fully optimized dynamics. Observe that in the unoptimized case, the population oscillations end somewhat short of the 
complete cycle, and the initially populated lcvel never achicves complete depopulation. Optimizing only the driving 
frequency does indeed maximize the population oscillations by inducing the complete depopulation of the initially 
populated level during oscillations, but at the same time the final population oscillation stops even further from the 
full cycle than in the unoptimized case. Finally, introducing the population transfer period correction alongside the 
driving frequency correction yields the required result: complete cycle of maximized population oscillations. 

2. Strong perturbation 

Increasing the driving perturbation intensity to somewhat greater value, erg = 0.25 (Fo = 6.11409 * 10~ 4 a.u..), 
the limitations of the analytical theory clcarly emerge. This is shown in Fig. 3: Fig. 3.a - Fig. 3.c respectively show 
the unoptimizcd, analytically optimized (according to Ea. H18l) and Ea. l|34|l ') and 'manually optimized' dynamics. The 
corresponding pulse duration times, aimed at producing three complete population oscillations, are T n — 847324 a.u., 
T opt = 901075 a.u. and T man = 884300 a.u.. 

Notice that in the analytically optimized case, Fig. 3.b, the initially populated level still fully depopulates, which 
indicates that even for this rather strong perturbation, the frequency correction Ea. (|18fl still stands strong. However, 
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2.a Unoptímized dynamics 
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2.b Se mi- optimized dynamics 
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2.c Fully optimized dynamics 
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FIG. 2: Significance of the total analytical correction (in driving frequency and total pulse duration) for the dynamics of a 
mildly disturbed system. In all plots, major oscillations are the populations of the two targeted levels (a and 0) whereas the 
minor oscillations are the population of the perturbing level (p). For the value of perturbation strength parameter ap p = 0.05, 
the loss of the final unoptimized population transfer amplitude amounts about 2%. Optimization reduces this loss to below 
0.05%. The difference between the pulse duration obtained from Standard 7r-pulse theory and the optimized value is 1.6%. 
Right hand-side plots present the details from the left hand-side plots. 



the corrected period, although closer to the correct value than in the unoptimized case, is still somewhat removed 
from the correct value. Unfortunately, this 'optimization error' cannot be remedied analytically. Remember that the 
analytical result Ea . (|34(l is obtained using only the first order approximation ÍEa. H^lJl ) to the full dynamical equations 
Ea. (|29|l . With perturbation as strong as in this case, the dynamical impact of the neglected elements of that equation 
begin to show. However, as shown in Fig. 3.c, the full cycle of oscillations can still be produced, but this additional 
correction to the pulse duration had to be found by hand, using the trial and error method. 



B. Population transfer 



The two final examples demonstrate the application of the developed optimization theory to the most interesting 
dynamical case regarding the coherent control: that of the single population transfer between the two targeted levels 
a and (3. As the validity and the limitations of the whole theory were already explored in the previous two examples, 
the following examples will only demonstrate the improvements to the population transfer that can be produced using 
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3.b Analytically optimized dynamics 
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FIG. 3: Limitations of the analytical optimization theory. Again, major oscillations are the population of the two targeted levels 
whereas the minor oscillations are the population of the perturbing level. For the value of perturbation strength parameter 
is now <Xa p = 0.25, which is just beyond the limiting value for the full applicability of the presented optimization procedure. 
Although the analytical correction improves the final population transfer from 92% to 97% (with pulse duration correction of 
6%), the theory presented in this paper cannot account for an additional 1.6% correction in the duration of the pulse that 
further increases the final population transfer to over 99.99%. 

the above results. 



1. Legitimate perturbation 

Again as in the previous section, the first example (Fig. 4) presents the dynamics of the system subjected to the 
external drive of limiting intensity. In this case, the perturbation strength parameter amounts cr^ p = 0.1, corresponding 
to the F = 4.07606 * 10~ 4 a.u.. Calculated population transfer times are = 211831 a.u. and T opt = 218483 a.u.. 

The unoptimized (dotted line) and the optimized dynamics (sòlid line) are plotted on the same graph to facilitate 
the comparison between the two. Only the dynamics of the two target levels is shown - the plot of the perturbing 
level's (p) dynamics is omittcd for the sake of clarity of the overall graph. Although the loss of the population transfer 
in the unoptimized case is not great, it neverthcless is noticeable. On the othcr hand, introducing the corrections for 
pulse frequency and pulse duration clearly improves the population transfer bringing it very close to 100%. 
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4.a Unoptímized dynamics 




FIG. 4: Applicability of the analytical optimization procedure to the maximization of the population transfer. Yet again, major 
oscillations are the population of the two targeted levels whereas the minor oscillations are the population of the perturbing 
level. For this limiting value of perturbation strength parameter of a^ p = 0.10, the optimization almost completely eradicates 
the loss of the population transfer due to the dynamical impact of the perturbing level, increasing the population transfer from 
96% to 99.7%. Optimized pulse lasts 3% longer than the one obtained from the Standard 7r-pulse theory. 



2. Extreme perturbation 

The final example - presented in Fig. 5. - is qualitative, rather than quantitative, but even as such it is quite 
indicative of the ovcrall uscfulncss of the whole optimization theory. The perturbation is now extreme, with strength 
parameter a^ p = 1 corresponding to the Fo — 1.22282 * 10 -3 a.u.. Calculated population transfer times are T v = 
70610 a.u. and T opt = 82816 a.u.. 

Again, the unoptimized and the optimized dynamics are plotted on the same graph. Although optimization now 
clearly does not lead to the complete population transfer, the improvement from the unoptimized dynamics is signif- 
icant demonstrating that even for this perturbation intensity the developed optimization theory qualitatively works 
quite nicely. 



11 



n(t) 


5.a Unoptímized d v na inics 








0.8 








0.6 


\ 






0.4 
0.2 


















t T 


n(t) 


5.b Optimized dynamics 
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FIG. 5: Breakdown of the quantitative optimization, but qualitatively the theory is still applicable. Half way through the pulse, 
in the unoptímized case now the population of the perturbing level surpasses the population of the initially unpopulated level. 
Introduction of the optimized driving frequency and pulse duration significantly - although not fully - rectifies the population 
transfer from below 40% to almost 90%. Optimized pulse is now almost 20% longer then the one obtained from the ordinary 
7r-pulse theory. 

IV. CONCLUSION 

The aim of research that led to this paper was to explore the possibility of using 'old fashioned' and rather simple 
phcnomenon of Rabi oscillations for the controllcd manipulation of the population in general many level system. This 
paper rounds up the tòpic of analytical optimization of pulse paramctcrs (frequency chirp and pulse duration) , opened 
in the author's previous work (Q) that would lead to maximizing the population transfer between two targeted levels 
of the system. The theory developed provides the exact quantitative predictions of to what extent the dynamical 
impact of the remainder of the many level system (beyond the two levels selected for the population transfer) begins 
to interfere with the targeted population transfer. It also provides the closed (albeit recursive) analytical expressions 
for the optimization of pulse parameters. 

Although the major correction to the population transfer is achieved by optimizing the driving pulse's frequency 
chirp (given in |9j ) , this paper provides the additional fine tuning by establishing similarly simple analytical expression 
for the determination of the optimal pulse duration. It demonstrates that the Standard formula of the 7r-pulse theory, 
Eq. (|19J) begins to fail as the perturbation increases to and beyond the well defined limiting value. It also provides 
some remedy to this failure. 



12 



The whole theory presented in [9( and this paper deals with only single làser pulse driving one particular transition 
in the many level system. The further research currently under way considers the possibility of applying a number of 
distinct but simultaneous optimized pulses to drive the population through the chain of transitions through the system, 
hence producing as clean as possible transfer between the two levels not coupled by the single photon transition. The 
prcliminary results indicate that an analytical optimization formula can be developed even for such a case. 



[1] G.K.Paramonov, V.A.Savva, Resonance effects in molecule vibrational excitation by picosecond làser pulses, Phys. Lett. 
A 97A (8) (1983) 340-342. 

[2] S.Chelkowski, A.Bandrauk, P.B.Corkum, Efficient molecular dissociation by chirped ultrashort infrared làser pulse, Phys. 

Rev. Lett. 65 (19) (1990) 2355-2358. 
[3] M. Kaluza, J.T.Muckerman, P.Gross, H.Rabitz, Optimaly controlled five-laser infrared multiphoton dissociation of hf, J. 

Chem. Phys. 100 (6) (1993) 4211-4228. 
[4] K.Bergmann, H.Theuer, B.W.Shore, Coherent population transfer among quantum states of àtoms and molècules, Rev. 

Mod. Phys. 70 (1998) 1003-1026. 
[5] H.Rabitz, Shaped làser pulses as reagents, Science (299) (2003) 525-527. 

[6] M.S.Shariar, P.Pradhan, Fundamental limitation on qubit operations due to the bloch-siegert oscillation, Proceedings of 
the quantum communication, measurement and computing (QCMC'02) quant-ph/0212121 

[7] J. CA. Barata, W.F.Wreszinski, Strong coupling theory of two level àtoms in periòdic fields, Phys. Rev. Lett. 84 (10) (2000) 
2112-2115, |quaïït -ph/9906029 

[8] K.Fujii, Two-level system and some approximate solutions in the strong coupling regime, quant-ph/0301145 . 

[9] D.Bonacci, S.D.Bosanac, N.Doslic, Analytic pulse design for selective population transfer in many-level quantum systems: 
maximizing the amplitude of population oscillations, Phys. Rev. A 70 (4) (2003) 043413, quant-ph/0310034 
[10] M.Holthaus, B.Just, Generalized tt pulses, Phys. Rev. A 49 (3) (1994) 1950-1960. 

[11] R. Nagle, E. B.Saff, Fundamentals of differential equations and boundary value problems, 2nd Ed., Addison Wesley, 1996. 



